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1.0  EXECUTIVE  SUMMARY 


Ocean  Systems  Research,  Inc.  under  the  sponsorship  of  the  Naval  Air  Development  Center  and 
the  Office  of  Naval  Technology  has  recently  completed  testing  modifications  to  their  Rapid  Thermal 
Ice  Penetrator.  The  tests  demonstrated  that  the  use  of  higher  energy  propellants  results  in  faster 
penetration,  reduced  propellant  volume/weight,  and  a  simplified  design. 

These  tests,  which  were  conducted  through  10  foot  thick  columns  of  ice,  clearly  substantiated 
the  performance  advantages  to  be  gained  by  increasing  the  energy  content  of  the  propellant.  Using 
the  same  amount  of  propellant  (21.3  pounds)  as  was  used  in  our  successful  Arctic  penetration  test,  the 
10  foot  ice  column  was  penetrated  in  approximately  100  seconds.  Since  the  total  bum  time  was  141 
seconds,  this  indicates  that  30  percent  of  the  propellant  can  be  eliminated,  while  still  achieving 
complete  penetration  through  10  feet  of  ice. 

We  predict  that  with  additional  increases  in  propellant  energy  content,  and  further  refinements 
to  the  nozzle  to  increase  melting  efficiency,  we  can  reduce  propellant  volume  by  more  than  40  percent 
as  compared  to  our  Arctic  tests. 

The  chart  below  highlights  chronologically,  the  technological  advancements  made  in  the  Rapid 
Thermal  Ice  Penetration  program  to  date,  along  with  a  comparison  with  the  competing 
thermochemical  approach. 

Rapid  Thermal  Ice  Penetrator 
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Although  this  contract  dealt  only  with  improvements  to  the  ice  penetrator  itself,  it  should  be 
noted  that  many  significant  technical  hurdles  have  already  been  overcome  in  the  development  of  the 
uprighting  device.  A  necessary  component  of  the  autonomous  ice  penetrating  sensor  package,  the 
uprighter  raises  the  entire  penetrator  from  horizontal  to  vertical  once  it  has  come  to  rest  on  the  ice, 
and  it  provides  guidance  for  the  penetrator  body  during  initial  penetration  into  the  ice.  The  prototype 
"A"  size  uprighter  devices  fabricated  for  Arctic  testing,  used  a  combination  of  high  tension  coil  and 
leaf  springs,  and  6  legs  that  encircled  the  penetrator  body,  with  a  maxinuiTn  thickness  of  .43  inches 
and  a  weight  of  less  than  10  pounds.  Successful  operation  of  the  device  was  demonstrated  (with  a  50 
pound  system  weight)  at  APLIS  91. 

The  uprighter  device  also  incorporated  a  revolutionary  remote  release  device.  No  existing 
release  mechanism  could  be  found  that  could  hold  the  200  pound  spring  force,  release  it  upon 
demand,  and  fit  within  our  size  constraints.  OSR  therefore  designed,  built,  and  tested  a  special 
release  device  that  weighed  less  than  an  ounce,  was  capable  of  holding  over  500  pounds,  yet  could  be 
instantaneously  released  by  a  small  9.6  volt  battery  pack.  The  device,  which  was  totally  inert  (i.e., 
contained  no  explosives  or  pyrotechnic  devices),  performed  reliably  in  every  Arctic  test.  This  remote 
release  technology  has  enormous  potential  for  any  application  requiring  extremely  high  release  forces. 
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2.0  BACKGROUND 


At  the  APLIS  91  ice  camp,  OSR,  under  the  sponsorship  of  the  Naval  Air  Development  Center 
and  the  Office  of  Naval  Technology,  successfully  demonstrated  the  ability  to  autonomously  penetrate 
10  feet  of  ice  in  under  two  minutes.  This  represented  a  major  milestone  towards  the  ultimate  goal 
of  this  program,  which  is  to  develop  an  air  deployable  ice  penetrating  sensor  package  that  fits  within 
the  standard  "A"  size  configuration  (i.e.,  36  inches  long  by  4  7/8  O.D.).  The  units  tested  in  the  Arctic 
included  a  fully  self-contained  uprighting  device  (see  Figures  lA  and  IB)  that  uprighted  the 
penetrator  from  horizontal  to  vertical  and  provided  initial  guidance  to  ensure  vertical  penetration. 

In  order  to  stay  within  the  size  constraints  of  the  "A"  size  configuration,  however,  and  still 
allow  a  sufficient  volume  for  the  sensor  payload,  the  volume  of  propellent  must  be  reduced.  Based 
upon  the  size  of  the  units  tested  at  APLIS  91,  this  reduction  in  propellant  volume  must  be  on  the 
order  of  30-40  percent.  The  most  straightforward  method  of  achieving  this  reduction  in  propellant 
volume  is  by  using  a  propellant  with  a  higher  energy  density.  Although  our  analytical  model 
indicated  that  increasing  the  energy  content  of  the  propellant  allows  a  reduction  in  propellant  volume 
for  the  same  depth  of  penetration,  this  had  never  been  tested.  All  of  the  previous  ice  penetration 
testing  performed  by  OSR  had  used  the  same  basic  propellant.  This  testing  effort,  therefore,  was 
undertaken  to  determine  if,  in  fact,  use  of  higher  energy  propellants  would  lead  to  the  significant 
reductions  in  propellant  volume  needed  to  meet  the  program  requirements. 

3.0  DESCRIPTION  OF  HARDWARE 

The  purpose  of  this  testing  effort  was  to  prove  that  increasing  the  energy  density  of  the 
propellant  would  singularly  lead  to  decreased  propellant  volume  requirements.  Therefore,  in  order 
to  keep  costs  low  and  to  limit  the  introduction  of  new  variables,  the  existing  design  and  (where 
possible)  hardware  was  used  for  this  testing.  Accordingly,  two  of  the  motor  cases  that  had  been  used 
in  the  Arctic  tests  were  refurbished.  In  addition,  two  nozzle  closures  used  in  the  Arctic  tests  were 
also  refurbished  by  burning  out  the  nozzle  inserts  and  insulator  material,  and  by  welding  over  the 
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Figure  1A  Uprighter  Device  in  Closed  Position 

(Shown  on  pallet  with  ice  penetrator  inserted,  prior  to 
installation  of  release  mechanism  band  clamp.) 


Figure  1B  Uprighter  Device  In  Open  Position 


existing  nozzle  openings.  Thus  the  external  envelope,  size,  and  configuration  of  these  test  units  was 
identical  to  those  tested  in  the  Arctic. 

To  reduce  the  potential  for  clogging  of  the  nozzles  (as  was  experienced  during  the  Arctic 
testing),  several  changes  were  made  to  the  internal  materials  and  configuration.  Post  test  analysis  of 
the  Arctic  units  indicated  that  glass  was  the  primary  constituent  of  the  substance  clogging  the  nozzles. 
Since  glass  was  a  major  component  of  the  insulator  used  for  the  Arctic  tests,  an  alternate  insulator 
material  was  selected  for  these  tests.  The  insulator  selected  was  a  Kevlar  based  EPDM  with  no 
appreciable  glass  content. 

The  internal  nozzle  configuration  and  materials  were  also  changed.  The  erosion  characteristics 
of  the  higher  energy  propellant  are  somewhat  different  than  the  propellants  used  in  earlier  testing. 
The  Arctic  units  utilized  a  molybdenum  nozzle  insert  whereas  the  new  units  used  a  4-D  carbon- 
carbon  matrix  insert.  This  change  was  required  for  several  reasons.  The  higher  energy  propellant 
bums  at  a  higher  flame  temperature,  which  is  above  the  operating  region  for  the  molybdenum.  In 
addition,  the  increased  particulate  content  of  the  higher  energy  propellant  exhaust  tends  to  erode  the 
molybdenum.  Lastly,  the  increase  in  percentage  of  aluminum  (fuel)  decreases  the  free  oxygen  in  the 
combustion  products  thereby  reducing  the  reactivity  of  those  gases  with  carbon.  One  added  benefit 
from  the  switch  to  the  carbon  insert  is  that  it  is  much  easier  to  machine/fabricate. 

Following  the  August  30  test,  inspection  revealed  significant  erosion  occurred  in  the  silica 
phenolic  insulator  used  in  the  nozzle  closure.  To  prevent  the  recurrence  of  this  erosion,  the  silica 
phenolic  insulator  was  burned  out  and  replaced  with  a  carbon  phenolic  for  the  second  test  (tested  16 
September). 

The  configuration  of  the  nozzle  insert  was  also  modified  to  eliminate  any  obstructions  in  the 
exhaust  flowpath  that  could  potentially  lead  to  clogging  of  the  nozzles.  Figures  2A  &  B  highlight 
these  differences.  Note  that  the  Arctic  units  (Figure  2A)  had  several  steps/ledges  where  material 
could  coagulate.  These  steps/ledges  were  eliminated  in  the  two  latest  unit  tested.  (See  Figure  2B.) 
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Figure  2A  Nozzle  Configuration  -  Arctic  Units 


Figure  2B  Nozzle  Configuration  -  High  Energy  Units 


3.1  Propellant  Selection 

The  average  energy  content  of  the  propellants  tested  prior  to  this  effort  was  approximately 
2300  Btu/Lbm.  Our  goal  for  this  testing  effort  was  to  select  a  propellant  with  an  energy  content  of 
over  3000  Btu/Lbm,  however,  due  to  other  design  and  performance  considerations  a  slightly  lower 
energy  propellant  was  selected.  The  propellant  chosen  for  these  tests  had  an  energy  content  of 
approximately  2900  Btu/Lbm. 

The  following  rationale  was  used  in  the  propellant  selection  process.  Figure  3  shows  several 
possible  propellant  formulations  suitable  for  this  application.  From  Figure  3  it  is  readily  apparent  that 
flame  temperature,  as  well  as  energy  content  (Btu/Lbm),  increases  almost  linearly  with  increasing 
aluminum  content  in  the  propellant.  Based  upon  this  chart  alone  it  would  seem  that  propellant  TP- 
3340  with  18  percent  aluminum  is  a  good  choice.  Among  the  other  factors  that  had  to  be  considered, 
however,  are  propellant  mass  flow  rate  and  chamber  pressure.  All  of  our  previous  tests  were  designed 
to  produce  a  mass  flow  rate  of  between  0.14  and  0.16  Lbm/sec  (see  Figure  4).  This  mass  flow  rate 
results  in  a  approximate  2  to  2  1/2  minute  burn  time  for  the  given  propellant  load.  From  Figure  4 
one  can  see  that  in  order  to  achieve  a  mass  flow  in  the  desired  operational  region  using  the  TP-3340 
propellant,  a  chamber  pressure  in  excess  of  1000  psia  would  be  required. 

Since  the  existing  cases  were  designed  for  a  500  psia  chamber  pressure,  a  1000  psia  chamber 
pressure  would  require  that  new  cases  be  fabricated.  The  higher  pressure  would  also  require  a  smaller 
nozzle  diameter  which  is  contrary  to  our  desire  to  reduce  the  potential  for  nozzle  clogging. 
Alternately,  operating  with  the  TP-3340  propellant  at  under  500  psia  results  in  a  mass  flow  rate  of 
approximately  0.10  Lbm/sec,  which  translates  into  a  burn  time  of  approximately  3.5  minutes  for  the 
21.0  lb  propellant  load.  This  burn  time  was  deemed  to  be  too  far  outside  the  2.0  minute  requirement. 
Propellant  TP-3396  was  therefore  chosen  because  its  mass  flow  rate  (at  300  psia)  provided  us  with 
the  desired  burn  time,  and  its  energy  content  of  approximately  2900  Btu/Lbm  was  high  enough  above 
the  previously  tested  propellants  to  provide  us  with  appropriate  performance  characteristics. 
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Figure  3  Propellant  Characteristics 
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Figure  4  Propellant  Mass  Flow  Comparison 


4.0  TESTING 


4.1  Test  #1 

The  first  of  the  two  ice  penetrator  tests  using  a  higher  energy  propellant  was  conducted  on 
30  August.  The  motor  contained  approximately  21.3  pounds  of  propellant  with  an  energy  content  of 
about  2900  Btu/Lbm.  It  had  one  forward  facing  0.313  diameter  central  nozzle,  and  two  0.85  diameter 
nozzles  diametrically  opposed  at  120*  for  reverse  thrust.  The  test  setup  was  similar  to  previous  tests 
(see  Figure  5)  in  that  the  penetrator  was  mounted  atop  a  3  foot  by  3  foot  by  10  foot  tall  ice  block. 
An  absolute  encoder  was  attached  to  the  rear  of  the  penetrator  body  for  penetration  distance  vs.  titnft 
measurements. 

Following  ignition,  the  motor  penetrated  approximately  2  feet  of  ice  in  just  over  20  seconds, 
then  appeared  to  stop  moving.  The  motor  went  on  to  bum  for  approximately  3  minutes  and  10 
seconds,  achieving  a  final  penetration  depth  of  about  4  feet.  Post  test  inspection  revealed  that  the 
forward  nozzle,  as  well  as  one  of  the  reverse  facing  nozzles,  had  become  totally  clogged.  The  one 
remaining  nozzle  opened  up  significantly,  thus  explaining  the  extended  bum  time.  The  material 
blocking  the  two  nozzles  was  found  to  be  aluminum  (a  constituent  added  to  the  propellant  to  increase 
flame  temperature).  A  detailed  review  of  the  video  tape  recording  of  the  test  indicated  that  clogging 
of  the  forward  nozzle  began  at  about  5  seconds  and  it  appeared  fully  clogged  around  20  seconds. 

4.2  Test  #2 

The  second  ice  penetrator  test  using  the  high  energy  propellant  was  conducted  on  16 
September.  The  motor  contained  21.32  pounds  of  propellant  with  an  energy  content  of  about  2900 
Btu/Lbm.  It  had  a  single,  forward  facing  central  nozzle  with  a  diameter  of  0.325  inches.  It  was 
designed  to  operate  at  a  chamber  pressure  of  282  psia  with  an  average  mass  flow  rate  of  0.152 
Lbm/sec.  Total  predicted  burn  time  was  140  seconds.  Approximately  30  pounds  was  added  to  the 
rear  of  the  penetrator  giving  it  a  total  weight  of  about  70  pounds.  To  prevent  the  motor  from  starting 
with  the  nozzle  flush  against  the  ice  surface,  an  extender  cup  was  welded  to  the  front  of  the  nozzle 
(see  Figure  6). 
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Figure  5  In-House  Test  Arrangement 

(Ice  Penetrator  mounted  in  hollow  tube  atop  10  foot  tall 
by  3  foot  square  ice  block.) 


6  Penetrator  Nozzle  With  Extender  Cup 


Prior  to  the  test,  two  cracks  were  noted  in  the  front  surface  of  the  ice,  but  their  cause  could 
not  be  determined  as  to  the  cause  of  the  cracks.  Two  possible  explanations  are  that  they  occurred 
during  transport  (via  forklift)  from  the  cold  chamber  to  the  test  site,  or  that  they  were  caused  by 
thermal  stresses  due  to  the  90+  degree  ambient  temperature.  Following  ignition,  the  motor  began 
penetrating  the  ice  immediately.  The  average  penetration  rate  was  approximately  1.2  inches  per 
second.  Penetration  occurred  at  about  1  minute  40  seconds,  and  the  motor  continued  to  bum  for 
approximately  2  minutes,  21  seconds.  The  resultant  hole  was  6  to  6.5  inches  in  diameter  and  appeared 
perfectly  vertical  (see  Figure  7). 

5.0  ANAL YnCAL  MODELING 

The  analytical  model  developed  by  OSR  under  a  previous  NADC  contract  was  used  to  help 
evaluate  the  various  nozzle  configurations  prior  to  fabrication  and  testing.  In  addition,  following  each 
test,  the  actual  results  were  used  to  improve  the  model’s  correlation. 

As  stated  previously  in  this  report,  the  penetrator  used  in  test  #1  had  one  central  downward 
facing  nozzle  and  two  reverse  acting  nozzles.  The  predicted  bum  time  was  131  seconds,  and  the 
predicted  thrust  was  approximately  23  pounds.  The  model  predicted  penetration  in  approximately 
95  seconds  (see  Figure  8A)  unfortunately,  the  nozzle  clogged  during  the  test  and  successful 
penetration  was  not  achieved. 

Test  #2  utilized  the  same  propellant  as  test  #1,  however  it  had  only  one  central  nozzle.  The 
predicted  thrust  was  approximately  28  pounds.  The  model  predicted  penetration  in  approximately 
110  seconds  (see  Figure  8B).  The  total  system  weight  of  the  Test  #2  unit  was  approximately  70 
pounds.  To  depict  the  effect  of  thrust  on  penetration  rate.  Figure  8C  shows  that  due  to  its  higher 
thrust,  the  Test  #2  unit  would  not  penetrate  10  feet  of  ice  if  its  weight  were  reduced  to  50  pounds. 

A  copy  of  the  computer  code  including  the  required  model  inputs  is  provided  in  Appendix 
A.  The  theoretical  ice  penetrator  performance  data  for  the  Test  #2  unit  is  provided  in  Appendix  B. 
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Figure  7  Test  #2  Results 

(Nearly  vertical  through  the  10  foot  thickness) 
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6.0  CONCLUSIONS/LESSONS  LEARNED 

Based  upon  this  successful  testing  using  a  higher  energy  propellant,  several  conclusions  were 
reached,  i.e.. 

Use  of  the  higher  energy  propellant  makes  an  "A"  size  sensor  package  very  feasible. 

A  single  central  nozde  eliminates  the  nozzle  clogging  problem. 

Spinning  is  not  required  for  vertical  penetration  over  the  10  foot  ice  thickness. 
The  average  penetration  rate  of  1.2  inches  per  second  is  our  fastest  to  date.  The  single  central  nozzle 
also  appears  to  be  the  most  efficient  design.  It  produced  the  smallest  diameter  (approximately  6 
inches),  and  the  smoothest  and  straightest  hole  to  date. 

Although  the  test  program  was  a  tremendous  success  there  are  still  some  areas  requiring 
additional  work.  Since  this  was  a  propellant  improvement  program  only,  no  improvements  to  the 
autonomous  uprighting  device  were  made.  With  regard  to  the  ice  penetrator  motors,  however,  we  feel 
very  strongly  that  virtually  every  problem/risk  area  has  been  worked  out,  save  one.  The  only 
remaining  effort  standing  in  the  way  of  fully  self  contained  "A"  size  ice  penetrator/sensor  package 
is  a  nozzle  study/test  program  to  reduce  the  thrust  from  over  20  pounds  to  approximately  5  pounds. 

There  are  several  options  available  for  reducing  the  thrust  including  adjusting  the  propellant 
properties,  varying  the  chamber  pressure,  adding  a  diffuser  to  the  nozzle,  etc.  It  is  anticipated  that 
the  final  low  thrust  design  will  incorporate  some  or  all  the  these  modifications.  Based  upon  the 
outcome  of  these  two  tests,  however,  it  is  clear  that  increasing  the  energy  density  of  the  propellant 
increases  the  amount  of  ice  melted  per  pound  of  propellant.  The  positive  effect  of  this  finding  is  that 
we  can  now  use  energy  release  rate  (i.e.,  Btu/Lbm)  rather  than  mass  flow  rate  as  a  primary 
characteristic  for  propellant  selection.  Since  Ais  no  longer  locks  us  in  to  the  0.14  to  16  Lbm/sec 
flow  rate,  even  higher  energy,  lower  bum  rate  propellants  now  become  excellent  candidates.  Of 
course,  the  lower  bum  rates  mean  lower  mass  flow  rates  which  translates  into  reduced  thrust.  And 
fortunately,  since  the  future  propellant  loads  will  be  shorter,  a  reduced  burn  rate  will  be  needed  to 
achieve  the  total  bum  time  required  (i.e.,  2-3  minutes). 
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APPENDIX  A 


Ocean  Systems  Research  Ice  Penetration  Hodel  (09-29-91  )(K.  Tzou) 
user's  information 


Input  data 

Ice(in)  a  Total  Depth  of  Ice 

U_tota((Lbs)  =  Initial  total  weight  of  the  system 

U_propel(Lbs)=  Initial  total  weight  of  the  propellant 

B_rate( Lbs/sec)  =  Burning  rate  of  the  propellant 

M_cjet  =  No.  of  center  jet 

M_dJ1  =  No.  of  down-ward  jet  type#1 

N_dj2  =  No.  of  down-ward  jet  type#2 

N_ujet  =  No.  of  up- ward  jet 

d_cjet(in)=  diameter  of  center  jet 

d_dj1  (in)=  diameter  of  down- ward  jet  type#l 

d_dj2  (in)=  diameter  of  down-ward  jet  type#2 

d_ujet(in)=  diameter  of  up-ward  jet 

ang_cjet(deg)=  angle  of  center  jet 

ang_dj1  <deg)=  angle  of  down-ward  jet  type#1 

ang_dj2  (deg)s  angle  of  down- ward  jet  type#2 

ang_ujet(deg)=  angle  of  up- ward  jet 

Dia<in)  3  diameter  of  the  ice  penetrator 

Height  (in)  =  height  of  the  ice  penetrator 

son_vel(fn/sec)=  sonic  velocity  (from  Thiokol) 

Mach  3  Mach  number  (from  Thiokol) 

Rho(g/cc)  3  exit  gas  density  (from  Thiokol) 
Vis(lbf-sec/ft**2)  3  exit  gas  viscosity  (from  Thiokol) 

Spe.  Heat(cal/gm-degK)3  exit  gas  specif  heat  (from  Thiokol) 

Pr  3  Exit  gas  mixture  Prandtl  nijiter(from  Thiokol) 

c^time,  c_gwo,  c_cro,  t_start  :  correlation  coefficients 


Output  Results 

time  (sec)  :  time  in  seconds 
Depth(in)  :  depth  in  inches 
Weight(Lbs):  total  remaining  weight 
v( in/sec)  :  rate  of  penetration 


1)  prepare 

an  input  file  called 

”sep_70.inp“. 

TP-H-3396  Propellant 

-  September  16,  1991 

_ * 

Ice  (in) 

-120.0 

► _ 

W^total 

70.0 

_ *_ 

W_Propel 

21.32 

Berate 

0.152 

N_cjet 

1.0 

N_dj1 

0.0 

_ _ 

M_dj2 

0.0 

N_uiet 

0.0 

d_cjet 

0.465 

d^djl 

0.000 

d_dj2 

0.000 

d^ujet 

0.00 

ang^cjet 

0.0 

f _  _ 

_ 

ang^djl 

0.0 

ang_dj2 

0.0 

ang^ujet 

0.0 

Oia(in) 

4.0 

_ 

Height 

44.9 

* 

son^vel 

934.6 

Mach 

2.048 

Rho  Vis 

5.0989e-4  0.1662e-5 

Spe. Heat 

0.48884 

Pr 

0.55906 

'2).  prepare  a  run  file  called  ”  go_ice-bat  ** 

copy  sep_70.inp  ice.inp 
osr^ice 

opy  ice. out  sep_70-out 
rase  ice. out 
print  sep_70.out 


[  5).  to  run  osr^ice  program  type 
go_i  ce 


illed  **  sep_70.out  **  will  be  printed  out. 

TP-H-3396  Propellant  -  September  16,  1991 

r  *  .  *  *  _ 

Ice  (in) 

-120.0 

U^total  W^Propel 

70.0  21.32 

_  _ *_ 

B_rate 

.152 

_ _ 

N^cjet 

1. 

M_dj1  H_dj2 

0.  0. 

_ * _ 

N_ujet 

0. 

.... _ ★ _ 

... _ ★ 

_ 

d^cj'et 

.4650 

_ _  *_ 

d.dji  d_dj2 

.0000  .0000 

d^ujet 

.0000 

ang^cjet 

0. 

ang_dj1  ang_dj2 

0.  0. 

_  _ _ *_ 

ang^ujet 

0. 

.  __ 

Dia(in) 

4.0 

Height  son^vel 

44.9  934.60 

...... • ..... . ..*. 

Hach 

2.048 

_ 

Rho 

.3099E-03  . 

Vis  Spe.Heat 

1662E-05  .4888 

Pr 

.5591 

....... 

C^time 

100.00 

c_gwo  c^cro 

.800  4.000 

*  * 

t^start 

5.000 

_ _ 

_ * 

.  _ _ 

time(sec) 

Depth(in)  Weight (Lb)  V( in/sec) 

5.48 

-1.00 

70. 

2.0675 

5.97 

-2.00 

70. 

2.0521 

6.46 

-3.00 

70.  ^ 

2.0368 

6.96 

-4.00 

70. 

2.0215 

7.46 

-5.00 

70. 

2.0061 

7.96 

-6.00 

70. 

1.9908 

8.46 

-7.00 

70. 

1.9755 

8.97 

-8.00 

69. 

1.9602 

9.49 

-9.00 

69. 

1.9449 

10.01 

-10.00 

69. 

1 .9296 

10.53 

-11.00 

69. 

1.9143 

11.06 

-12.00 

69. 

1.8991 

11.59 

-13.00 

69. 

1.8838 

12.12 

-14.00 

69. 

1.8686 

12.66 

-15.00 

69. 

1.8533 

13.20 

-16.00 

69. 

1.8381 

Ocean  Systems  Research  Ice  Penetration  Model  (09-29-91)(IC.  Tzou) 
Program  osr^ice 

real  k_gas,N_jet,N_ujet,N_^djl ,N_dj2,N_cjet,Mu_ave,k_ice, ice_Z,Mach 
character  title*60,  tit1*60,  tit2*60 
dimension  Z_totaU150),  T_total<150),  t_sec<150) 
dimensiorv  h_gas<150),  h_water(150),  Ueight(150) 

OPEN  (UMIT=s1,  FILE=M'ce•inp^  STATUS='OLD') 

OPEN  (UNIT=2,  FILE=M'ce-OUt',  STATUS='NEWO 


Input  data 

lce(in)  ss  Total  Depth  of  Ice 

U_total(Lbs)  =  Initial  total  weight  of  the  system 

U j3ropel(Lbs)=  Initial  total  weight  of  the  propellant 

B_rate(lbs/sec)  =  Burning  rate  of  the  propellant 

M_cjet  »  No.  of  center  jet 

M_dJ1  No.  of  down- ward  jet  type#1 

N_dj2  s  No.  of  down- ward  jet  type#2 

N_ujet  =  No.  of  up- ward  jet 

d_cjet(in)«  diameter  of  center  jet 

d_dj1  (in)=  diameter  of  down- ward  jet  type#1 

d_dj2  (in)=  diameter  of  down-ward  jet  type#2 

d_ujet(in)s  diameter  of  up- ward  jet 

ang_cjet(deg)=  angle  of  center  jet 

ang^djl  (deg)=  angle  of  down- ward  jet  type#1 

ang_dj2  (deg)^  angle  of  down-ward  jet  type#2 

ang_ujet(deg)=  angle  of  up-ward  jet 

Oia(in)  =  diameter  of  the  ice  penetrator 

Height  (in)  *  height  of  the  ice  penetrator 

son_veUm/sec)s  sonic  velocity  (from  Thiokol) 

Mach  =  Mach  number  (from  Thiokol) 

Rho(g/cc)  =  exit  gas  density  (from  Thiokol) 
Vis(lbf-sec/ft**2)  =  exit  gas  viscosity  (from  Thiokol) 

Spe.  Heat(cal/gm-degK)=  exit  gas  specif  heat  (from  Thiokol) 

Pr  =  Exit  gas  mixture  Prandtl  number (from  Thiokol) 

c_time,  c_gwo,  c_cro,  t^start  :  correlation  coefficients 


Output  Results 

time  (sec)  :  time  in  seconds 
Depth(in)  :  depth  in  inches 
Ueight(Lbs):  total  remaining  weight 
v( in/sec)  :  rate  of  penetration 


read  input  data 


read  (1,100)  title 
write(2,200)  title 
print  100,  title 

100  format (a60) 

200  format(/20x,a60/) 

201  format (20x,a60) 

read  (1,102)tit1,tit2, 
write(2,202)tit1,tit2, 
read  (1 ,102)tit1,tit2, 
write(2,203)tit1,tit2, 
read  (1,102)tit1 ,tit2, 
write(2,209)tit1,tit2, 
read  (1,102)tit1,tit2, 
write(2,203)tit1,tit2, 
read  (1 , 102)tit1,tit2, 
write(2,204)tit1,tit2. 


ice_Z,  W_total, 
ice_Z,  W^total, 
N_cjet,  N_dj1, 
N_cjet,  N_dj1, 
d_cjet,  d_dj1, 
d_cjet,  d^djl, 
ang_cjo,ang_dj1 
ang_cjo,ang_dj1 
Dia,  Height, 
Dia,  Height, 


W^propel,  Berate 
W jjropel,  B_rate 
Njdj2,  M_ujet 
N_dj2,  N^ujet 
d_dj2,  d_ujet 
d_dj2,  d_ujet 
,  ang^dj  2 , ang_u j  et 
,  ang_d j  2 , ang_u j  et 
Son^vel,  Mach 
Son^vel,  Mach 


»  read  (1J06)tit1,tit2,  Rho,  Vis,  Cp,  Pr 

write<2,206)tit1,tit2,  Rho,  Vis,  Cp,  Pr 

I  read  (1,102)titl,tit2,  C_tiine,  c_gwo,  c_cro,  t^start 
I  write(2,208)tit1  ,tit2,  C^time,  c_gwo,  c^cro,  t^start 
read  (1,100)tit1 
r'  write(2,201)tit1 
I  102  format(a60/a60/5f10.0) 

^  202  fonnat(20x,a60/20x,a60/20x,2f10.1,  f10-2,  f10.3) 

^  203  fonnat(20x,a60/20x,a60/20x,4(f8.0,2x)) 

I  204  fonnat(20x,a60/20x,a60/20x,2fl0,1,f10,2,  f10.3) 

I  106  format(a60/a60/2e10.4,2f10.0) 

206  fonnat(20x,a60/20x,a60/20x,2e10.4,2f10.4) 

208  fonnat(20x,a60/20x,a60/20x,f10.2,3f10.3) 

209  format(20x,a60/20x,d60/20x,4f10.4) 

210  format(20x,a60/20x,a60/20x,2f10.4,3fl0.2) 
write(2,220) 

220  f onnat(/20x, '  time<sec)  Depth(in)  Weight(Lb)  V(in/sec)') 

*  -  assumption  &  initial  conditions 

2_totaU1)=0.0 

T^tota I ( 1 )=t_start 

t_sec(1)  =0.0 

Cm=0.85 

pi=3. 1415927 

d_c  j  et=d_c  j  et/ 12.0 

d_dj1  =d_djl  /12.0 

d^dj2  =d_dj2  /12.0 

d_u  j  et=d_u  j  et/ 12.0 

a_c  j  et=(pi  /4 .  >’*d_c  j  et**2 

a_d  j  1  =<  pi /4 . )*d_d j 1 **2 

a_dj2  =(pi/4,)*d_dj2*^2 

a_ujet=(pi/4.  )*d_ujet**2 

aN_c  j  et =N_c  j  et*a_c  j  et 

aN^djl  =M_djl  *a_dj1 

aN_dj2  =N_dj2  *a_dj2 

aN_ui  et=N_u  j  et*a_u  j  et 

aN_  j  et=aN_c  j  et+aN_d  j  l+aM_d  j2+aN_u  j  et 

c_d1=(90.  -ang^dj  1  )^c_gwo/90.0 

c_d2=(90.  -ang_d  j2)*c_gwo/90.0 

C_gw=  ( aN_c  j  et  *c_g  wo+aN_d  j  1  *c_d1 +aN_d  j  2*c_d2 )  /aN_  j  et 
:kt  print  *,  'C_gw  *',C_gw 

cr_dl =(90 . 0-  ang_d  j  1  )*c_cro/90 . 0 
cr_d2=(90 . 0  -  ang_d  j  2  )*c_cro/90 . 0 

C_crack=(aN_cjet*c_cro+aN_dj1*cr_d1+aN_dj2*cr_d2)/aN_iet 
:kt  print  *,  'C_crack=' ,C_crack 

C_Nu=0.268 
C_Re=0.625 
C_u=6.630 
k_ice=1.25 
alpha=0.0450 
T=32.0 
Tf=212.0 
Ti=-15.0 

cu_ang=cos<  pi  *ang_u  j  et/180 .0) 
cdl_ang=cos(pi *ang_di 1/180.0) 
cd2_ang=cos (pi *ang_d j  2/180.0) 

Rho_steam=1 .86 
Dia=0ia/12.0 

*  - - - - assume  initial  diameter  of  ice  hole  is  3”  bigger 

*  than  the  diameter  of  the  penetrator 
D_i  ce=0  i  a-**3 .0/12.0 


i  360  formate/'  Time  (sec)  =' ,f6.1,5x, 'Depth  (in)  s'^f6.0) 
close  (1) 
close  (2) 
stop 
end 


Following  files  are  in  this  disk 


Volune  in  drive  B  has  no  label 
Directory  of  B:\ 


P  FILE 

OAT 

598 

9-29-91 

12:39p  — 

----> 

this  file 

I  SEP_70 

OUT  ‘ 

10581 

9-29-91 

11:07a  — 

----> 

sainple  output  file 

*  SEP_70 

INP 

1048 

9-29-91 

10:52a  — 

sample  input  file 

OSR^ICE 

FOR 

7621 

9-29-91 

12:27p  — 

fortran  program 

1  DSR^ICE 

EXE 

36912 

9-29-91 

11:06a  — 

- > 

conpiled  run  file 

1  ^EAD 

ME 

13682 

9-29-91 

12:24p  — 

- > 

user's  information 

GOJCE 

BAT 

101 

9-29-91 

12:27p  — 

batch  run  file 

7  File(s)  1141760  bytes  free 


APPENDIX  B 


THEORETICAL  ROCKET  PERFORMANCE  ASSUMING  EQUILIBRIUM  COMPOSITION  DURING  EXPANSION 
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I  Height=Height/12-0 

* . -  assume  crack  area  *  C^crack  times  circular  area 

r  area^ice  =(pi/4.)*  D_ice**2 

[  area^cro  =C_crack*  area^ice 

areajsen  =(pi/4.)*  pia**2 
p  .  Weight(1)^_total 

B_t  i  me=W jDrope  I  /Berate 

M  J  et=M_c  j  et+N_d  j  1 +N_di  2 

s_  j  et=sqrt  ( pi  *D  i  a**2/  ( 4 .  *N_  j  et ) ) 

I  u__ex  i  t =Cm*Son__ve  l*Mach*3.2808 

Lxt  print  */u_exit\  u^exit. 

k_gas=Cp*Vi  s*32 . 2*3600 .  /Pr 
!ct  print  */k_gas^k_gas 

z_del=1. 0/12.0 
Rho=Rho*62.43 

. assune  50%  gas  50%  water  at  212  deg  F 

Rho^mi x=0 . 5*Rho+0 . 5*Rho_steam*32 . 2 
ckt  print  */Rho_mix' ,Rho_mix 

B_force=0.0 

. - — . — . do 

do  300  k=2,140 

Z_total(k)=Z_total(k-1)  -  z^del 
Z_inch=Z_tota I < k)*1 2. 0 

if(-Z  totaUk).LE. Height)  B_force=-Z_totaUk)*areaj>en*Rho_mix 
if (-Z_total(k).GT. Height)  B_force=Height*areaj3en*Rho_mix 
kt  if  (k.EQ.2)  print  */B_force',  Bjorce 

I  Wei gh t ( k ) =Wei gh t <  k- 1 )  - t_sec( k- 1 )*B_rate 

T_Wei  ght=Wei  ght  ( k )  -  B_f  orce 

F_uiet=N_u  jet*Rho*a_u  jet*u_exi  t*C_u*u_exi  t*d_u  jet*cu_ang 
F  dji  =N  dji  *Rho*a_dj1  *u_exit*C_u*u_exit*d_dj1*cd1_ang 
F  dj2  =N  dj2  *Rho*a_di2  *u_exit*C_u*u_exit*d_dj2*cd2_ang 
F_c  jet=N_c  jet*Rho*a_c  jet*u_exi  t*C_u*u_exi  t*d_c  jet 
zJet=<F_c]et-F_uiet+F_d]1+F_d]2)/T_weight 
Re_ave=C_u*  ( u__ex  i  t  *d_c  j  et /z_  J  et  )*s_  j  et*Rho/ V  i  s 
Nu_ave=C_Nu*Re_ave**C_Re 
h_gas(  k  )=(  C_Nu*k_gas/s_J  et  )*Re_ave**C_Re 


kt  if  (k.EQ.2)  print  */h_gas',h_gas(k) 

* _ assume  h_water=C_gw*h_gas 


h_wa t er ( k ) =C_gw*h_gas ( k ) 

- - - -assume  h_water  vary  with  time- .C_time 

t i me_sec=T_tota  U  k- 1 ) 

h_water(k)=h_water(k)*exp(-time_sec/C_tinie) 
ratio^h  =h_water(k)/h_water(2) 


area^cra  «  area^cro  *  1.0 
c11-Log((T-Tf)/(Ti-Tf)) 

time=-kjce*z_del*areajce*c11/(alpha*h_water(k)*  are&^^ra) 
t_sec ( k) =t i me*3600 . 0 
Ve  I  j3en*2_de  I  /  t_,sec  ( k ) 

Ve  l_i  ps=Ve  l_pen*  12.0 
T_totaUk)=T_totaUk“1  )+t_sec(k) 
zin=z_jet*12.0 

kt  if(k.EQ.2)  print  */2jet(in)=',zin 
if<T_totaUk).GT.B_time)  go  to  320 
if  (Zjnch.LT.icej:)  go  to  320 
write(2,260)T_total(k),2_inch,Weight(k),VelJps 

261  format(2f10.2) 

260  format(20x,f10.2,4x,f10.2,4x,f8.0,7x,f8.4) 

300  continue 

320  print  360  ,  T.totaUk-1),  ZJnch+1.0 


